Constraints on the IR behavior of the ghost propagator in Yang-Mills theories 
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We present rigorous upper and lower bounds for the momentum-space ghost propagator G{p) of 
Yang-Mills theories in terms of the smallest nonzero eigenvalue (and of the corresponding eigenvec- 
tor) of the Faddeev-Popov matrix. We apply our analysis to data from simulations of SU(2) lattice 
gauge theory in Landau gauge, using the largest lattice sizes to date. Our results suggest that, in 
three and in four space-time dimensions, the Landau-gauge ghost propagator is not enhanced as 
compared to its tree-level behavior. This is also seen in plots and fits of the ghost dressing function. 
In the two-dimensional case, on the other hand, we find that G{p) diverges as p~^~^'^ with k « 0.15 , 
in agreement with Ref. [3| . We note that our discussion is general, although we make an application 
only to pure gauge theory in Landau gauge. Most of our simulations have been performed on the 
IBM supercomputer at the University of Sao Paulo. 

PACS numbers: 11.15.Ha 12.38.Aw 14.70Dj 



I. INTRODUCTION 

In the Gribov-Zwanziger scenario 0], confinement of 
quarks (in Landau and in Coulomb gauge) is related to 
a ghost propagator G{p) enhanced in the infrared (IR) 
limit when compared to the tree- level propagator 
An enhanced ghost propagator is also required by the 
Kugo-Ojima criterion Q in order to explain confinement 
of color charge (in Landau gauge). 

In order to investigate the origin of this IR enhance- 
ment we consider, following Ref. [4|, a generic gauge con- 
dition T[A\ = 0, where A is the gauge field. The con- 
dition is imposed on the lattice by minimizing a func- 
tional E[U], where U is the (gauge) link variable. From 
the second variation of E[U] one can define the Faddeev- 
Popov (FP) matrix M{a,x;b,y). (Here we use a and 
b to indicate color indices, whereas x and y represent 
points of the lattice.) In the SU(7Vc) case, A4 is an 
{N^ - 1)V X {N^ - \)V matrix, where ]/ = iV^ is the 
lattice volume, given by the number of points iV along 
each side of the lattice and by the space-time dimension 
d. This matrix is real and symmetric with respect to 
the exchange (a, a;) (&, y) (see ReL @ for a thorough 
discussion of the properties of the FP matrix in Landau 
gauge). At any (local) minimum of £■[[/], all the eigenval- 
ues of M. are positive (modulo trivial null eigenvalues). 
The set of all minima of i?[J7] is the so-called Gribov 
region 17. On the boundary of — known as the first 
Gribov horizon 9f7 — the smallest nontrivial eigenvalue 
Amin of the FP matrix is null. Since the configuration 
space has very large dimensionality we expect that, in 



the limit of large lattice volumes V , entropy favors con- 
figurations near i.e. \min should go to zero 
in this limit. This is indeed the case in 2d d, 3d i and 
4d Landau gauge in 4d Coulomb gauge [§[ and in 
4d Maximally Abehan gauge (MAG) 11). 

The ghost propagator is written in terms of the inverse 
of the FP matrix as 



G(p) 



-It:! k-{x~y) 



V 



■M ^{a,x;a,y), (1) 
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where the lattice momentum p{k) has components 
Pf,{k) = 2 sin {iTkf,/N) with fc^ = 0, 1, . . . , TV - 1. Also, 
here and in the formulae below, a path-integral average 
is understood. Since M{a,x;b,y) develops a null eigen- 
value at the Gribov horizon dfl , one might expect that 
the corresponding ghost propagator diverges at small mo- 
menta in the infinite- volume limit. This could in turn 
introduce a long-range effect in the theory, related to the 
color-confinement mechanism Q. Indeed, in several 
numerical studies (using relatively small lattice volumes), 
the ghost propagator G{p) seems to be IR-enhanced in 
3d and 4d Landau gauge [H, and in 4d Coulomb 
gauge [ij]. On the other hand, in MAG [ll| one finds an 
IR-finite G{p). At the same time, recent results in Lan- 
dau gauge using very large lattice volumes [l^ [l^ 113 
suggest an essentially tree-level-like IR ghost propagator 
G{p) in 3d and in 4d. (A similar result has been obtained 
in Refs. [H, using different analytic approaches.) Fi- 
nally, in the 2d case, one finds [l| G{p) ^ p~^'^ after 
extrapolating data for the IR exponent k, to infinite vol- 
ume. Therefore, we have cases for which G{p) is not 
IR-enhanced in the infinite- volume limit and thus the ar- 
gument reported above cannot be valid in general d- 
The aim of this work is to understand under what condi- 
tions one should expect an IR-enhanced ghost propagator 
G{p). This is done by investigating constraints on the be- 
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havior of G{p) and thus obtaining better control over its 
infinite-volume extrapolation. Clearly, since the only di- 
verging quantity (besides the lattice volume) is 1/Xmin, 
the IR enhancement of G{p) may be closely related to 
this eigenvalue and to the projection of its corresponding 
eigenvector ^min{a,x) on the plane waves e"^'^*'^ '^/'^!/ 
for small momenta p{k). 

II. LOWER AND UPPER BOUNDS FOR G{p) 



TABLE I: The ghost propagator G(ps) for the smallest 
nonzero momentum ps = 2sin(7r/A'^) and the inverse of the 
smallest nonzero eigenvalue Xmin of the FP matrix for var- 
ious lattice volumes and P values in the 3d case. Data (in 
physical units) are taken from Ref. [H. 



N 


P 


a-' [GeV] 


G{p,) [GeV-"] 


1/A™„, [GeV-"] 


20 


6.0 


1.733 


7.3(1) 


32.5(7) 


30 


6.0 


1.733 


19.4(3) 


94(2) 


20 


4.2 


1.136 


21.0(4) 


107(3) 


30 


4.2 


1.136 


54.5(8) 


282(6) 



In this section we obtain upper and lower bounds for 
the momentum-space ghost propagator G{p) in terms of 
the smallest nonzero eigenvalue Xmin [and of the cor- 
responding eigenvector ^minia, x)] of the FP matrix 
M{a,x;b,y). To this end, let us first introduce our no- 
tation. As usual, we define the generic eigenvalue A; of 
A4{a,x;b,y) and the corresponding eigenvector ipi^a^x) 
by using the relation Ai{a, x;b,y) ipi{b,y) = Xit/ji^a^x), 
where the index i takes values 1,2,..., {N^ ~ 1)V and 
the sum over repeated indices is understood. As stressed 
in the Introduction, this matrix is positive (or semi- 
positive) definite, whenever we use a minimizing condi- 
tion in order to fix a gauge (minimal gauge). Also, since 
it is obtained from a second-order expansion, this matrix 
can always be written in a symmetric form. Then, the 
eigenvectors '0i(a, x) can be assumed orthogonal to each 
other and normalized as ^ 2:) ^/^j (a, x)* — Sij, 
where * indicates complex conjugation. Using this nota- 
tion we can write 



M{a,x;b,y) 



E 



■4't{a,x) ipiib, y)* 



(2) 



Note that we are working in the space orthogonal to the 
kernel of the FP matrix. Then we can write the inverse 
of this matrix as [lO] 

M-\a,x-b,v) = E A-iV.(a,x)V.(6,y)* • (3) 
By using Eq. JT]) we obtain 

G^p) = E E^r'i^^.Kp)P, (4) 

where tpi{a,p) — '^^x^ii^^^) I is the Fourier 

transform of the eigenvector ^i{a,x). Since all the 
nonzero eigenvalues are positive and < Xmin < A^, we 
can write the inequalities 

1 1 



I'x E < G(P) , (5) 



i,\i^O a 



(Note that we assume nondegenerate eigenvalues. How- 
ever, the second inequality applies also when Xmin is de- 
generate and the first one can be easily modified to take a 



degeneracy into account.) By summing and subtracting 
in Eq. ([6]) the contributions from the eigenvectors corre- 
sponding to a null eigenvalue and using the completeness 
relation i^iia^x) ipi{b,y)* = Sab^xy, we find 



Gip) < 



1 



m - 1 



E Ei^^(«'^') 



(7) 

Let us stress that the above results are simply a conse- 
quence of the positivity of the FP matrix, i.e. they apply 
to gauge-fixed configurations that belong to the interior 
of the first Gribov region. 



III. BOUNDS IN LANDAU GAUGE 

In the Landau case the eigenvectors corresponding to 
the zero eigenvalue are constant modes, i.e. they con- 
tribute only to the case p = 0. Thus, for any nonzero 
momentum we have 



1 



1 



m - 1 A. 



J2 \^rmn{a,p)\^ < G{p) < 



1 



A. 



(8) 



As said in the Introduction, in the infinite-volume limit 
the measure gets concentrated on the boundary of the 
(first) Gribov region, i.e. Xmin goes to zero in this limit. 
In particular, as for the Laplacian operator, one can ex- 
pect to find Xmin ^ A^-" for large N. Now, if we make 
the hypothesis of a power-law behavior p"2-'2K j^j, G{p) 
in the IR limit and consider the smallest nonzero mo- 
mentum on the lattice [i.e. Ps — 2 sin(7r/A'^)], we have 
G{ps) ~ JSf^+^i^ in the limit of large N. Then from Eq. 
([5]) we get 2 + 2k < a. Thus, assuming a power-law 
behavior for G{p), a > 2 is a necessary condition to ob- 
tain a ghost propagator G{p) enhanced in the IR limit 
compared to the tree- level behavior 1/p^. 

The lower bound in Eq. ^ clearly depends on the 

behavior of the quantity \ipmin{a-,Ps)\'^ ■ Note that this 
quantity cannot diverge in the infinite- volume limit. (Ac- 
tually, since the eigenvectors are normalized, this quan- 
tity is always between and 1.) Thus, if we furthermore 
make the assumption of a behavior at large N given by 
N'"', with 7 > 0, we find the condition a — 7 < 2 -I- 2k. 

As a check of the proposed bounds, we collect here 
results for the ghost propagator G{ps) at the smallest 
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FIG. 1: The ghost propagator G{pa) for the smallest nonzero 
momentum ps — 2sin(7r/A'^) (in GeV~^) as a function of the 
inverse lattice side 1/L (GeV) for the 2d case (top, at /3 = 10 
with volumes up to 320'^), 3d case (center, at /3 = 3, with 
volumes up to 320"^) and the 4d case (bottom, at /3 = 2.2 
with volumes up to 128^). Data are taken from Ref. [l^ for 
the 3d and the 4d cases. We also show (in 2d and in 4d) 
the inverse of the smallest nonzero eigenvalue Xmin of the 
FP matrix (in GeV^^). In these two cases one can verify 
the inequality 1/Xmin > G{ps)- The fitting parameters are 
reported in Table Hill Note that we did not evaluate Xmin for 
all our configurations. 



TABLE II: As in Table U for the 4d case. Data (in lattice 
units) are taken from Refs. _^,'l2'] using the average called first 
copy. We find similar results when considering data obtained 
using the average called absolute minimum. 



N 


P 


G{ps) 


1 / Amin 


8 


0.8 


8.94(8) 


91(3) 


12 


0.8 


22.1(6) 


210(20) 


16 


0.8 


40.6(7) 


360(60) 


8 


1.6 


7.2(1) 


61(3) 


16 


1.6 


32.1(3) 


220(20) 


8 


2.7 


3.4(1) 


11.0(4) 


12 


2.7 


7.1(3) 


25(1) 


16 


2.7 


12.9(6) 


45(4) 



TABLE III: Fits of G{ps) and of 1/Xmin respectively using the 
Ansatze bL^^^"' and cL". In the 4d case the fit for 1/Xmm 
has been done by considering only the three largest physical 
volumes. When considering the three smallest physical vol- 
umes we find a = 0.9(3). In the other cases all data have 
been considered for the fit. 



d 


b 


2 + 2k 


c 


a 


2 


0.026(1) 


2.251(9) 


0.12(3) 


2.20(4) 


3 


0.11(3) 


2.02(4) 






4 


0.086(3) 


2.043(8) 


1.2(1) 


1.53(2) 



nonzero momentum = 2sin(7r/iV) and for the smallest 
nonzero eigenvalue Xmin from Refs. 0, IH, UBl for the 
3d and 4d cases. These data are shown for the smaller 
lattices in Tables U [IT] and for the larger lattices in Fig. 
[T] We also present new data in the two-dimensional case. 
(These data were obtained together with the data for the 
gluon propagator reported in AH data refer to the 

SU{2) case. Let us recall that recent studies [13, HI] have 
verified the analytic prediction that Landau-gauge gluon 
and ghost propagators in SU(2) and in SU(3) are rather 
similar. Thus, we expect that the analysis presented here 
should apply also to the SU(3) case. In all cases the 
quantity \/\min has been evaluated using the average 
value for Xmin and propagation of errors. 

As one can see, the upper bound in Eq. ([5]) is always 
satisfied. (Of course, this is the case also when p ^ ps.) 
The bound has also been verified using the data at /3 = 
from Refs. We have fitted the data for G{ps) and 

for Xmin as a function of L for the data shown in Fig. [fus- 
ing a power-law Ansatz. It is interesting to note that (in 
2d and in 4d) we actually find a smaller than 2 + 2k. The 
bound is still satisfied, since the multiplicative constant 
in the fit of Xmin is larger. It would be interesting to see 
how the upper bound is realized when even larger lattice 
volumes are considered. Also, note that the upper bound 
in the 4d case seems to saturate as the volume increases. 
This would indicate that the contribution to G{p) from 
all the eigenvalues Xi > Xmin stays finite at large vol- 
ume and that the exponent 7, defined above, is zero. Let 
us note that this is not the case at smaller lattice vol- 
umes, since then one needs to consider the contribution 
to G{p) from the first 150-200 smaller eigenvalues in or- 
der to reproduce the behavior of the propagator at small 
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TABLE IV: Fits of p^Gip^) using the Ansatze ap"^" in the 
2d case and a — 61og(l + cp^) in the 3d and 4d cases. In the 
first case we used the data in the range £ [0, 0.5] for the 
fit. In the other two cases the data in the range G [0, 1] 
have been considered. 



d 


a 


K 


b 


c 


2 


1.134(7) 


0.177(2) 






3 


4.7(1) 




0.579(5) 


320(20) 


4 


4.28(1) 




0.69(2) 


33(3) 



TABLE V: Fits of p^G{p^) using the Ansatze a(p"^* + 
&/'=)/(l+p^^) in the 2d case and a-&[log(l + cp^) + dpV(l + 
p^) in the 3d and 4d cases. In the three cases we used the 
whole range of momenta for the fit. 



d 


a 




6 


e 


c 


d 


2 


1.24(3) 


0.16(2) 


0.86(3) 


0.75(15) 






3 


4.75(1) 




0.491(5) 




450(30) 


7.1(1) 


4 


4.32(2) 




0.38(1) 




80(10) 


8.2(3) 



momenta 10]. Finally, let us recall that finite-size effects 
for G{p) (at fixed p) are generally very small. This can 
be explained considering the spectral density of the FP 
matrix, as done in Ref. 4]. Thus, our results indicate 
that in 3d and in 4d the ghost propagator is likely not 
IR-enhanced, while in 2d we obtain k « 0.1. We note 
that our result in the 2d case is essentially in agreement 
with Ref. where however the value of k is obtained 
by first fitting G{p) in an appropriate momentum range 
for each lattice volume and then performing the extrap- 
olation to infinite volume by a four-parameter fit. The 
fits for Xmin and G{ps) as a function of L are usually 
simpler, as shown above. 

These results are confirmed if one considers the dress- 
ing function p^G{p^) at the largest lattice volumes for 
the three cases (see Fig. [5]). Indeed, the data in the 2d 
case can be fitted using a power-law Ansatz ap~'^'^ (see 
Table HvT) . with k = 0.177(2), essentially in agreement 
with pj] . The fitted value for k decreases if one considers 
fewer points, dropping points with larger p, and becomes 
0.152(7) when considering data with p^ e [0,0.1]. The 
same Ansatz does not work in 3d and in 4d. In these 
cases the data for the dressing function can be described 
by the Ansatz a — b log(l -I- cp'^), inspired by the logarith- 
mic corrections suggested in Ref. [19]. In both cases one 
finds (approximately) p^G{p^) — > 4.5 in the limit p ^ 0. 
This result supports the value k « obtained above, 
when considering the bound using Xmin- It is interest- 
ing that one can also obtain a relatively good fit for the 
dressing function in the whole range of momenta by con- 
sidering the fitting function a{p~'^^ + hp'^^)l{\ + p^'^) in 
the 2d case and a - 6[log(l + cp^) + dp'^]/[l + p^) in 3d 
and in 4d. These fits and the corresponding parameters 
are reported in Fig. [5] and in Table |Vl 
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0.001 0.01 0.1 1 10 100 

FIG. 2: The ghost dressing function p^G{p^) as a function 
of (in GeV) for the 2d case (top, at /3 = 10 with volume 
320^), 3d case (center, at f3 — 3, with volume 240'^) and the 4d 
case (bottom, at /3 = 2.2 with volume 80*). The fitting func- 
tions and the corresponding fitting parameters are reported 
in Table lYl 



IV. CONCLUSIONS 



Our data suggest that the Landau-gauge ghost propa- 
gator does not diverge faster than 1 /p^ at small momenta 
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in 3d and in 4d, while in 2d G{p) does not diverge faster 
than p"^"^" with k between 0.1 and 0.2, in agreement 
with Ref. [ij. These results have been obtained by con- 
sidering the inequality in Eq. ([7]) and by fitting the data 
for G{ps) and Xmin a-s a function of L. In particular, we 
note that the use of the upper bound to constrain G{p) 
is quite convenient, since Xmin does not depend on p. It 
might thus be of interest to optimize the evaluation of 
and of ipmin{a,x). 
The results obtained here in the 3d and 4d cases, to- 
gether with those reported in Ref. '2l'| for the gluon prop- 
agator, seem to contradict the Gribov-Zwanziger confine- 
ment scenario for Landau gauge. On the other hand, one 
should establish if the behavior of gluon and ghost prop- 
agators at very small momenta is really so essential for 
the explanation of confinement. After all, when dynam- 
ical quarks are considered [l^, string breaking should 
manifest itself at a scale of about a fermi, i.e. for a en- 
ergy scale of about 200 MeV. Thus, for the physics of 
hadrons the behavior of the propagators at intermedi- 
ate momenta is probably more important. Let us recall 
that for a space-time separation of about 1 fermi the 
transverse gluon propagator violates reflection positivity 
[2^, i.e. becomes negative in real-space coordinates, and 



the exponent k for the ghost propagator is about 0.2-0.3 
considering p ~ 0.5 GeV [1, [H, [ij, [H [M [13 • Thus, 
in the range 200-500 MeV, nonperturbative effects are 
clearly present in the behavior of the gluon and of the 
ghost propagator. For this range of momenta the use of 
our lower bound for the ghost propagator could probably 
also be valuable. Of course, the open problem is if these 
effects suffice to explain color confinement. 
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